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jj 要 : Owens Valley Radio Observatory (OVRO)40 米 望远镜 的 观测 数据 显示 了 BL Lac 

天 体 CGRaBS J0141-0928 在 15 GHz 射电 波段 的 剧烈 光 变 . 我们 利用 Lomb-Scargle- 

Periodogram (LSP) 方 法 和 加 权 小 波 Z 变换 (Weight Wavelet 7Z-transform，WW2) 方 法 以 及 

Jurkevich 方法 分 析 了 CGRaBS J0141-0928 的 光 变 周期 ， 结 果 显 示 该 天 体 具 有 约 为 649 天 的 准 

周期 振荡 (Quasi-Periodic 0Scillation，QP0)， 其 置信 度 为 4. 46， 其 准 周 期 振荡 可 能 是 由 喷 

流 的 螺旋 运动 产生 . 我 们 用 双 指 数 函 数 拟 合 了 其 中 2 个 爆发 的 过 程 ， 从 而 得 到 其 光 变 时 标 ， 

r 进而 估算 其 平均 多 普 勒 因子 为 3. 8， 这 表明 CGRaBS J0141-0928 在 射电 辐射 波段 存在 显著 的 聚 
= 束 效应 .采用 DCF(Discrete Correlation Function) 方法 对 射电 波段 分 别 与 Gamma-ray. 36 
村 学 R 波段 的 光 变 相关 性 进行 分 析 ， 发 现 光 学 光 变 与 射电 光 变 有 强 相 关 性 ， 且 光学 光 变 超前 于 


射电 光 变 66 + 40K. 
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N 0 引言 
c WEARS (Blazar) 是 活动 星系 核 (Active Galactic Nuclei, AGNs) 中 性 质 较为 特殊 的 子 
N 类 ， 其 相对 论 性 喷 流 几乎 正 对 着 地 球 。 耀 变 体 具 有 极端 的 观测 特性 ， 包 括 高 光度 、 高 偏振 、 快 


速 光 变 ， 以 及 从 射电 到 高 能 y 射线 的 非 热 连续 辐射 。 耀 变 体 的 两 个 子 类 分 别 是 平 谱 射 电 耀 

ARS (FSRQ) 和 蝎 虎 座 BL 型 天 体 (BL Lac) ”。 其 中 BL Lac 天 体 的 光谱 只 有 一 些微 弱 的 发 射 

eU 线 或 者 没有 发 射线 ， 但 是 有 很 强 的 X 射线 及 7 射线 辐射 ”。 研 究 发 现 ， 一 些 BL Lac 天 体 的 长 

一 期 尺度 变化 是 周期 性 的 ， 且 不 同 波段 之 间 可 能 存在 不 同 的 联系 。 通 过 对 BL Lac 天 体 光 变 的 观 

£ 测 和 研究 能 获得 天 体 的 物理 机 制 ， 辐 射 过 程 和 内 部 结构 参数 等 重要 信息 ”“， 对 活动 星系 核 的 探 
索 和 了 解 有 重要 意义 。 

CGRaBS J0141-0928 是 一 颗 红 移 为 0. 733 的 炮 变 体 ”， 分 析 BL Lac 天 体 光 变 周 期 的 方法 有 

很 多 ， 最 常用 的 有 自 相 关 函 数 分 析 及 周期 拟 合法 ， 时 间 序 列 的 功率 谱 分 析 方 法 ，Jurkevich X 

法 ， 这 些 方法 都 十 分 广泛 的 应 用 在 BL Lac 天 体 的 长 周期 光 变 观 测 与 研究 中 “"。 本 文 分 别 使 用 

Lomb-Scargle periodogram 方法 ， 加 权 小 波 Z 变换 方法 和 Jurkevich 方法 对 CGRaBS J0141- 

0928 天 体 15 GHz 射电 波段 的 光 变 周期 进行 研究 。 利 用 双 指 数 函 数 拟 合 了 爆发 过 程 ， 并 估算 了 

多 普 勒 因子 。 使 用 离散 相关 函数 (DCF) 方 法 分 析 了 y 和 射电 波段 ， 以 及 光学 和 射电 波段 的 相关 

HET, rh LSP 方法 和 加 权 小 波 Z 变换 法 是 初次 用 于 研究 CGRaBS J0141-0928 天 体 的 光 变 周期 。 
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这 三 种 方法 都 非常 适合 处 理 非 等 间隔 数据 ， 得 到 的 结果 是 比较 稳定 的 周期 。 


1 样本 和 光 变 曲线 

美国 欧文 斯 谷 射电 天 文 台 (Owens Valley Radio Observatory, OVRO)fJ 40m 望远镜 是 全 
球 观测 兆 变 体 理想 的 设施 (https://sites.astro. caltech. edu/ovroblazars/)。 如 图 1 是 来 自 
于 (OVR0) 40 米 望 远 镜 BL Lac 天 体 CGRaBS J0141-0928 在 15 GHz 射电 波段 的 光 变 曲线 ， 其 中 
有 575 个 数据 点 。 从 光 变 曲线 可 以 看 出 : CGRaBS J0141-0928 在 射电 波段 的 活动 非常 激烈 ， 流 
量 随 着 时 间 波 动 起 伏 ， 除 了 在 2010 年 、2013 年 、2015 年 、2017 年 和 2019 年 有 五 次 明显 的 大 
爆发 外 ， 还 可 以 看 到 有 一 些 不 同 程度 的 小 爆发 ( 光 变 曲线 的 尖峰 处 ) 。 我 们 通过 计算 光 变 幅度 
Amp” 来 判断 天 体 的 活跃 程度 : 


Amp =100x,|(A,., -Ann Y —256’ percent. un 


其 中 As 和 Am 分别 表示 流量 的 最 大 值 和 最 小 值 ，Amp 值 越 大 表示 天 体 越 活跃 。 计 算得 到 PKS 
0139-097 天 体 在 射电 波段 的 光 变 幅度 为 62. 7， 表 明 CGRaBS J0141-0928 是 个 非常 活跃 的 天 体 。 
于 受 观测 条 件 的 限制 ， 光 变 曲 线 数据 不 连续 ， 因 此 周期 性 分 析 是 受到 限制 的 。 
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Fig.1 The light curves of CGRaBS J0141-0928 at radio band 


2 周期 分 析 
2. 1 Lomb-Scargle periodogram (LSP) 方法 

Lomb-Scargle periodogram (LSP) 是 一 种 广泛 用 于 寻找 准 周期 振荡 的 方法 CLomb1976; 
Scargle 1982; Press et al.1992) , HH Lomb 发 展 ， 经 过 Scargle 进一步 改进 。LSP 方法 不 
仅 能 在 一 定 程 度 上 减弱 时 域 序列 的 不 均匀 性 产生 的 虚假 信号 ， 还 可 以 有 效 地 从 时 域 序列 中 提取 
出 弱 周 期 信号 ”。 因 此 ，LSP 方法 能 很 好 的 寻找 隐藏 在 噪声 中 的 准 周期 振荡 光 变 。LSP 方法 的 
基本 原理 是 基于 傅立叶 变换 ， 将 一 系列 三 角 函 数 的 线性 组 合 = a cos ot + b sin ot i 
过 最 小 二 乘法 来 拟 合 时 间 序 列 ， 将 信号 特征 从 时 域 转换 到 频 域 上 。 基 本 公式 如 下 "": 
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ps ; oem ii “=f " Ss ii -| 


No No 
2r olt, -r) 2. em ext, —t) (2) 


其 中 +t 为 对 应 时 间 上 的 相位 修正 ， 计 算 公 式 为 : 


tan(2@r) = + 


Y cos(2a@r, ) (3) 


LSP 方法 在 对 正弦 型 和 数据 间隔 差距 较 小 的 时 间 序 列 进行 周期 分 析 时 有 较 好 的 处 理 效果 。 
为 了 验证 LSP 方法 计算 准 周期 的 正确 性 ， 我 们 首先 对 周期 图 进行 容 率 拟 合 ， 得 到 容 率 指数 B 
如 图 2。 然 后 我 们 计算 了 CGRaBS J0141-0928 的 置信 和 度 ， 结 果 如 图 3 所 示 。 图 3 中 ， 绿 色 线 表 
示 准 周期 图 ， 其 峰值 为 准 周期 结果 ; 蓝 色 、 红 色 、 紫 色 线 分 别 代 表 蒙 特 卡 罗 模 拟 的 95%、99%、 
99. 7% 置 信 度 。 图 3 中 的 绿 线 有 一 个 明显 的 峰值 为 649 天 ， 并 且 峰 的 置信 和 度 超 过 了 99. 7%， 说 
明峰 值 结 果 是 可 靠 的 ， 所 以 取 649 RAJ 1.78 年 作为 CGRaBS J0141-0928 天 体 的 准 周期 结果 。 
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图 2 CGRaBS J0141-0928 FACES Fait BEIGE RE 


Fig.2 CGRaBS J0141-0928 object logpower spectrum in radio band 
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图 3 CGRaBS J0141-0928 天 体 在 射电 波段 的 蒙特 卡 罗 模 拟 分 析 结 果 


Fig.3 CGRaBS J0141-0928 object Monte Carlo simulation analysis results in radio band 


2.2 ”加权 小 波 Z 变换 方法 


小 波 分 析 (Wavelet Aanlysis) 方法 是 一 种 新 的 分 析 方 法 ， 是 纯粹 数学 和 应 用 数学 的 结合 。 
Morlet 小 波 是 高 斯 包 络 下 的 单 频 复 正弦 函数 ， 是 一 种 复 小 波 。 利 用 Morlet 小 波 作 为 小 波 母 函 


数 进行 变换 ， 根 据 复 Morlet Kh: 


vit) _ e (gv = 


其 中 wo。 是 衰减 因子 。 
当 w。 取 较 大 值 时 ，Morlet 小 波 简化 为 : 


w(t) — e 12 oit, 


经 过 伸缩 a 和 平移 0 变换，(5) 式 变 为 : 


a 


ext) = etm (t-b)-cw,? (1-b . 


其 中 Wi=Wo/a, c-1/2W. 


CS pfi 
«(t8)-« a] 


(4) 


(5) 


(6) 


(7) 


根据 向 量 投影 的 方法 ， 可 以 认为 (7) 式 是 一 种 加 权 映 射 ， 它 以 6 (t) = e ”为 基 函 数 ， 
Wa se C “为 统计 权重 因子 。 另 外 引入 一 个 常数 函数 L(t) = 1, 可 以 得 到 向 量 空间 的 3 


个 基 函 数 : 
$,(t) = L(t), (8) 
p(t) = cos(o,( t — b)), (9) 
b(t) = sin(o,( t — b)) , (10) 
将 数据 向 量 x (t) 投影 到 这 3 个 基 函 数 上 , 得 : 
(11) 


3 
y= > y, AQ) 


Ya => sip]. (12) 


其 中 sc al Pb. 


根据 上 述 过 程 ，Foster 定义 了 加 权 小 波 变换 (WWT) |: 
(Na a IV, 


WWIee—— (13) 
2V, 
Kho (epo XEes6€)[Xes)| Dort) [Eet] 
Na, = EV oes a V, = z 
e Zo, 2.0, 2.0, 2 2.0, 2,0, 


由 于 小 波形 状 产生 变化 ， 低 频 部 分 的 有 效 数 据 Nore 会 多 于 高 频 部 分 的 有 效 数 据 ， 因 此 WT 
的 值 会 偏向 高 频 部 分 ， 使 得 结果 出 现 偏差 。 于 是 Foster 根据 Z 统计 量 定义 了 加 权 小 波 Z 变换 
P5 


I (14) 


它 满足 F 分 布 ， 自 由 度 为 3 和 2， 期 望 值 为 1。 利 用 (14) 式 可 以 得 到 如 图 4 的 周期 性 图 ， 
在 WZ 图 中 ， 频 率 轴 上 的 WZ 值 表示 数据 向 量 的 周期 性 ， 时 间 轴 上 的 WZ 值 则 表示 数据 向 量 随 
时 间 的 波动 情况 ”。 根 据 WZ 图 的 极 大 值 可 以 确定 天 体 光 变 曲线 的 周期 。 图 中 的 蓝 色 、 红 色 、 
紫色 线 分 别 代表 95%、99%、99. 7% 置 信 度 。 从 图 4 中 可 以 得 到 CGRaBS 70141-0928. 天 体 射 电波 
段 的 光 变 周期 约 为 636 天 ， 并 且 置 信 度 超过 了 99. 7%。 
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图 4 WWZ 方法 对 CGRaBS J0141-0928 在 射电 波段 的 分 析 结 果 


Fig. 4 WZ test result for period search in OGRaBS J0141-0928 at radio band 

2.3 Jurkevich 方法 
Jurkevich 方法 “由 Jurkevich 提出 于 1971 年 ， 是 一 种 基于 期 待 值 的 均 方差 的 周期 算法 ， 
通过 测试 周期 来 折 盖 数据， 适用 于 观测 数据 不 均匀 的 天 体 光 变 周 期 。 假 设 有 WNW 个 观测 的 样本 数 
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据 ， 每 次 的 测量 值 为 钱 了 为 所 有 样本 的 平均 值 , “为 测量 数据 样本 的 方差 ，S “为 标准 差 ， 则 
有 : 
Eo le (15) 2_ Wy aes. (16) 2 V? aD 
X=—> X,, V3 X? -NX?, ae 
wx 2 Pew 


根据 测试 周期 附近 的 相位 把 数据 样本 划分 为 mn 组 ， 则 第 1 组 的 统计 参数 为 : 


_ 1 N m; " y? 
X -2—Wxy. 9 V-WVvxogasx O9 s?=—i_, (20) 
i m, 之 i i 25 i m; i? i m, 一 1 
对 应 于 m 组 的 总 方差 为 : 
Ws V V?, (21) 
" 


当 测 试 周期 接近 于 真实 周期 时 ，V ;得 到 极 小 值 。 
此 外 ，Kidger 等 人 "在 Jurkevich 方法 的 基础 上 给 出 了 判断 周期 可 靠 性 的 方法 ， 即 : 
jeg 22 

广 :是 归 一 化 之 后 的 值 。 通 党 在 f 三 0. 25 时 ， 表 明光 变 可 能 具有 周期 性 ， 当 Z 20.5 时 ， 
表明 有 非常 显著 的 周期 。 

通过 Jurkevich 方法 ， 对 天 体 PKS 0139-097 射电 波段 光 变 曲线 的 分 析 结 果 如 图 5。 图 5 
中 的 绿色 线 以 下 代表 £F 20.25 时 的 置信 和 度 ， 红 色 线 以 下 代表 f£ m0.5 时 的 置信 度 。 寻 找 周 期 
时 ， 一 般 满足 以 下 两 个 条 件 的 周期 比较 可 靠 ， 一 是 分 析 的 数据 样本 时 间 跨 度 大 于 周期 的 六 倍 ， 
二 是 作出 的 曲线 有 明显 的 振幅 。 从 图 5 中 可 以 看 出 ， 对 应 于 VV; 的 最 小 值 ，CGRaBS J0141- 
0928 的 射电 波段 可 能 存在 650 天 的 比较 可 靠 的 周期 ， 对 应 的 £ £873 0.37. Jurkevich 方法 获 
得 的 650 天 周期 的 结果 和 LSP 方法 以 及 WWZ 方法 的 结果 非常 接近 。650 天 以 后 的 极 小 值 的 周期 
结果 明显 不 满足 以 上 两 个 条 件 ， 所 以 不 予 采 纳 。 
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图 5 Jurkevich 方法 对 CGRaBS J0141-0928 在 射电 波段 的 分 析 结 果 


Fig.5 The normalized Jurkevich test result for period search in CGRaBS J0141-0928 at radio band 
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3 多 普 勒 因子 分 析 

我 们 从 射电 的 光 变 曲线 中 挑选 出 2 个 包含 上 升 和 下 降 阶 段 的 爆发 过 程 ， 分 别 为 MJD 时 间 
56167. 4-56410. 7 和 56762. 8-57228. 5。 我 们 使 用 双 指 数 函 数 爆发 拟 合 公 式 …” 来 拟 合 这 2 个 
爆发 过 程 。 双 指数 函数 公式 如 下 ; 


tg-t tly zl 
rO-s+H。 ý ve] : (23) 


公式 (23) 中 ， 玉 表示 基底 流量 ， 加 表示 峰 对 应 的 时 间 ， 女 和 妃 分 别 表 示 指 数 上 升 和 下 降 的 
时 标 ，F。 衡 量 爆发 的 幅度 。 

双 指 数 函 数 拟 合 的 2 个 爆发 过 程 如 图 6 所 示 。 每 个 爆发 过 程 的 拟 合 参数 如 表 1， 第 1 列 为 
约 化 儒 略 日 的 范围 ; 第 2 列 为 拟 合 所 得 的 约 化 最 小 残 差 平 方 和 ; 第 3 列 为 拟 合 峰值 对 应 的 约 
化 儒 略 日 及 误差 ;第 4 列 为 基底 流量 及 误差 ， 第 5 列 为 爆发 的 幅度 及 误差 ， 第 6 列 和 第 7 列 分 
别 为 指数 上 升 时 标 和 下 降 时 标 及 误差 ， 第 8 列 为 每 个 爆发 过 程 所 对 应 的 多 普 勒 因子 。 男 外 ， 在 
图 6 中 的 右 图 中 可 以 看 到 一 个 额外 的 小 峰 ， 这 很 可 能 是 由 于 相对 论 喷 流 中 的 激 波 产生 的 光斑 导 
Skin. 

多 普 勒 因 子 (8 ) 与 喷 流 中 物质 流速 度 和 视角 有 关 ， 但 是 这 两 个 量 都 是 不 可 直接 观测 的 量 ， 
故 用 间接 方法 来 估算 多 普 勒 因子 是 有 必要 的 ， 其 中 通过 射电 光 变 估算 多 普 勒 因子 (5% 相对 比较 
准确 ””。 假 定 光 变 化 是 内 豪 的 ， 基 于 光 变 时 标 所 限定 的 源 的 线 度 大 小 ， 文 [30] 给 出 了 炮 变 体 
亮 温 度 的 计算 公式 : 


二 aD | (24) 
T, = (4.5x10 jar| s K 
AP, LATHE, AF ADA Jy 计 的 流量 变化 值 ， 刀 是 以 天 计 的 光 变 时 标 ， 和 是 以 厘米 计 的 观 
测 波长 ，D 是 以 百 万 秒 差距 (Mpc) 计 的 光度 距离 (采用 宇宙 学 参数 Hy =72 km s Mpc Ail 
Qu=0. 3) 。 

一 般 认为 焰 变 体 的 喷 流 成 分 对 应 的 亮 温度 7 不 能 超过 平衡 亮 温 度 75x 10"K ^U, Beit 
算出 来 的 亮 温 度 为 772.574x10"K 和 7F2.923x10 尿 ， 超 过 平衡 亮 温度 几 个 量 级 ， 表 明 存 在 显 
著 的 多 普 勒 增 亮 效 应 。 类 似 文 [28] ， 我 们 选择 7 作为 内 豪 亮 温度 。 利 用 如 下 公式 ， 进 一 步 估 
算 了 射电 波段 的 光 变 多 普 勒 因子 : 


[至 | l (25) 
T. 


eq 


Iu 


H (25) 式 得 到 多 普 勒 因子 65 的 值 分 别 为 3. 72 和 3. 88， 平 均值 约 为 3.8, Fan 等 人 ”估算 出 了 
y-ray 的 多 普 勒 因子 56=5. 50， 其 结果 和 我 们 估算 的 多 普 勒 因子 比较 接近 。 
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图 6 双 指 数 函 数 对 CGRaBS J0141-0928 的 2 个 爆发 (flare) 过 程 的 曲线 拟 合 


Fig.6 The double exponential function fits the curves of the 2 burst processes of CGRaBS J0141-0928 object 


31 爆发 过 程 的 拟 合 结 及 多 普 勒 因子 
Table1 The results of flare fitting and the Doppler factor 


1 1 1 1 L 1 1 1 iJ 
56600 56700 56800 56900 57000 57100 57200 57300 57400 


peak B t. F. Fo t. ta On 
56167.4- ^ 5.55e-04 56315.94+20.26 0. 43-0. 03 0.9750. 14 52. 583-21. 08 76.74-27.60 3.72 
56410. 7 
56762.8- ^ 1.82e-03 57028. 54+12. 97 0. 39-0. 06 0. 72+0. 06 34. 927. 61 94.98424.74 3.88 
57228. 5 


4 相关 性 分 析 


时 引入 的 ， 这 种 方法 可 以 通过 计算 时 间 延 迟 来 研究 天 体 的 结构 和 其 他 性 质 ”。 
离散 相关 函数 方法 的 定义 为 :如 有 任意 两 个 离散 数据 的 序列 a 和 bj ， 则 离散 相关 函数 的 
值 为 


UDCF, = (a, —a)x (^, -b) 
O, XO, (26) 


式 中 ，a 和 b 分 别 为 数据 序列 a 和 bi 的 平均 值 ，c .和 o ,分 别 为 相应 的 标准 偏差 。 
每 一 个 数值 UDCF5 与 延 时 TY=Ar=trti 有 关 。 而 对 于 有 噪 的 数据 ， 我 们 可 以 使 用 式 子 


(oz eZ Nor 一 亿 ) 来 代替 上 式 中 的 ,和 cv。 对 一 给 定 的 ,如 果 有 AM 个 UDC iE T-A +/2 
< Atul THA 1/2, 然后 将 这 M 个 数据 点 求 平均 值得 到 : 


离散 相关 函数 分 析 法 (DCF) 可 以 用 来 分 析 两 组 离散 的 数据 的 相关 性 “”。 该 方法 不 需要 对 
数据 样本 做 任何 的 处 理 就 可 以 判断 出 两 组 数据 的 相关 性 。DCF 方法 是 由 Edelson 研究 时 间 延 迟 


202111.00042v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1 
DCF(r)= ap 22UDCR,. (27) 

DCF ( t ) 就 是 离散 相关 函数 ， 在 延 时 域内 对 延 时 进行 区 间 划 分 ， 每 一 区 间 间 隔 为 Ar ， 这 
样 我 们 就 可 以 得 到 一 个 非常 有 用 的 DCF (+ ) 。 如 果 在 某 个 区 间 没 有 数据 点 ， 则 DCF (t) BUB e 
当 两 个 相关 序列 不 相同 时 ， 可 以 得 到 离散 互相 关系 数 。 当 两 个 相关 序列 相同 时 ， 则 得 到 离散 自 
相关 系数 。 在 离散 相关 函数 的 分 析 图 中 ， 如 果 DCF 的 峰值 越 大 ， 则 说 明 两 列 数据 的 相关 性 越 强 ， 
反之 则 越 弱 。 当 DCP 的 峰值 在 大 于 零 的 一 侧 时 , 表示 数据 a 超前 于 数据 p， 当 DCF 的 峰值 在 小 
于 零 的 一 侧 时 ， 则 表示 数据 a 滞后 于 数据 b。 通 过 离散 相关 函数 对 CGRaBS J0141-0928 射电 波 
段 分 别 和 来 自 Fermi Gamma-ray 空间 望远镜 y 波段 、 来 自 KAIT Fermi AGN Light-Curve 
Reservoir (http://herculesii. astro. berkeley. edu/kait/agn/) 光学 了 波段 的 数据 进行 相关 
性 分 析 ， 运 用 Fortran 编写 程序 计算 结果 如 图 7、8， 图 中 的 顶部 是 射电 波段 和 y 波段 /光学 R 
波段 的 流量 图 ， 底 部 是 射电 波段 和 y 波段 /光学 R 波段 的 相关 性 结果 图 ， 峰 值 越 接 近 1 说 明 相 
关 性 越 好 。 所 以 从 图 7 中 可 以 看 到 y 和 射电 波段 的 DCF 最 大 值 为 0.3， 相 关 性 非常 弱 ， 这 表 
明 它 们 的 辐射 区 域 不 同 并 且 辐 射 过 程 也 是 不 一 致 的 ， 图 8 中 显示 光学 和 射电 波段 的 DCP 最 大 值 


达到 了 0.71， 光 学 波段 超前 于 射电 波段 16~110 天 ， 说 明 它们 的 辐射 过 程 是 一 致 的 。 
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图 7 DCF 方法 对 CGRaBS J0141-0928 在 射电 和 y 波段 的 相关 性 分 析 。 顶 图 中 红色 和 黑色 分 别 表示 Gamma-ray 和 射电 波段 的 光 变 曲线 ， 

左边 标注 表示 射电 波段 的 单位 ， 右 边 标注 表示 Gamma-ray 波段 的 单位 。 底 图 表示 Gamma-ray 和 射电 波段 之 间 的 相关 性 结果 即 DCF fü. 
Fig. 7 Correlation analysis of CGRaBS J0141-0928 in radio and y bands by DCF method 
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图 8 DCF 方法 对 CGRaBS J0141-0928 在 射电 和 光学 R 波段 的 相关 性 分 析 。 顶 图 中 蓝 色 和 黑色 分 别 表示 Optical 和 射电 波段 的 光 变 
线 ， 左 边 标 注 表示 射电 波段 的 单位 ， 右 边 标 注 表 示 Optical 波段 的 单位 。 底 部 表示 Optical 和 射电 波段 之 间 的 相关 性 结果 即 DCF 值 。 
Fig.8 Correlation analysis of CGRaBS J0141-0928 in radio and optical R bands by DCF method 


5 讨论 与 结论 


通过 收集 BL Lac 天 体 CGRaBS J0141-0928 的 光 变 数据 ， 运 用 Lomb-Scargle 
periodogram 方法 ， 加 权 小 波 Z 变换 法 以 及 Jurkevich 方法 对 其 射电 波段 的 光 变 曲线 的 周期 
进行 分 析 ， 分 析 的 结果 相互 印证 。Lomb-Scargle periodogram 方法 分 析 射 电波 段 得 到 的 周 
期 约 为 649 天 ， 加 权 小 波 Z 变换 法 分 析 射 电波 段 得 到 的 周期 大 致 为 636 K, Jurkevich 方法 得 
到 的 周期 大 致 为 650 天 。 加 权 小 波 2 变换 法 和 Jurkevich 方法 的 计算 结果 进一步 支持 了 
CGRaBS J0141-0928 天 体 的 射电 波段 有 一 个 可 靠 的 约 为 649 天 左右 的 光 变 周期 。 利 用 双 指 数 函 
数 我 们 拟 合 了 其 中 2 个 爆发 的 过 程 ， 估 算 了 多 普 勒 因子 p= 3.8。 结 果 表 明 CGRaBS JO141- 
0928 天 体 射 电波 段 的 辐射 存在 显著 的 聚 束 效 应 ， 支 持 了 相对 论 性 喷 流 模型 。 通 过 离散 相关 分 
析 法 对 射电 波段 和 Gamma-ray、 射 电波 段 和 光学 R 波段 分 别 进行 相关 性 分 析 ， 结 果 显 示 Gamma- 
ray 和 射电 波段 之 间 存 在 弱 的 相关 性 ， 光 学 R 波段 和 射电 波段 之 间 存 在 比较 强 的 相关 性 ， 并 且 
光学 R 波段 超前 于 射电 波段 16 110 天 。 对 于 类 星体 长 周期 光 变 的 物理 模型 现在 仍然 不 是 很 清 

， 为 了 解决 这 些 问 题 提出 了 一 些 物 理 模 型 。 常 见 的 物理 模型 有 双 黑 洞 模型 ””" 、 螺 旋 喷 流 模 
=" ge ROR Ba 性 等 。 

CGRaBS J0141-09287 天 体 产 生 的 准 周期 可 能 是 由 喷 流 的 螺旋 运动 ”产生 的 。 喷 流 的 螺旋 
进 动 由 超大 质量 黑洞 〈SMBBH) 系统 的 轨道 运动 驱动 。Rieger 给 出 了 实际 物理 驱动 周期 PRU 
测 准 周期 P 之 间 的 关系 公式 为 : 


P 
z (28) 


1+ 
其 中 yy, 是 体 洛 伦 兹 因子 ， 约 等 于 7.5”，z 是 红 移 。 我 们 使 用 649 天 的 准 周期 得 到 实际 物理 驱 
动 周期 Pe. 71yr。 另 外 ， 如 果 使 用 y=15 作为 参数 ”， 可 以 得 到 P2230. 85yr。 主 黑洞 和 
次 黑洞 之 间 的 质量 比 为 〖 三 1/3 时 ， 称 为 “Major merger" ( 主 合 并 ) 。 若 质 量 比 为 3x 


P, Rm 


AN di 


minaX IV 合 


10'"， 则 称 之 为 “Minor merger”【〔 次 合并 )“”。 无 论 质 量 比 为 多 少 ， 主 黑洞 的 质量 都 可 以 ) 
以 下 公式 估计 


M ~ PPHOY0 Mo (29) 

其 中 以 年 为 单位 ， 对 于 SMBBH 系统 的 主 并 合 ， 质 量 比 可 以 假设 为 〖 =3/2。 将 参数 代入 (29) 
SO ee E el dr e 
并 合 的 SMBBH 系统 的 主 黑洞 质量 为 M x10””Mo。Wu 等 人 给 出 了 CGRaBS J0141-0928 天 体 的 


黑洞 质量 为 MLO Me， 其 结果 和 我 们 估算 的 主 黑洞 质量 一 致 。 
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Analysis of optical property of BL Lac object CGRaBS 
J0141-0928 
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Relativistic Astrophysics, Nanning 530004; 3 China; School of Physical Science and Engineering Technology, Guangxi University, Nanning 
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Abstract: Observation data from the Owens Valley Radio Observatory (OVRO) 40-meter 
telescope showed multiple dramatic light variations in the 15 GHz radio band of the blazar CGRaBS 
J0141-0928. LSP (Lomb-Scargle-periodogram) method, WWZ method (Weighted Wavelet Z- 
transform, WWZ) and Jurkevich method are used to analyze the light variable period of CGRaBS 
J0141-0928. The results show that the object has a light variable period of about 649 days (4.46 
significance level), and its quasi-periodic oscillation (QPO) may be caused by the helical motion of 
the jet. The brightness temperature, time scale of light variability and Doppler factor were estimated 
by fitting the processes of two bursts with double exponential function. The Doppler factor estimated 
by the optical time scale of variability is 3.48, and the results show that CGRaBS J0141-0928 has a 
significant beaming effect in the radio radiation band. The DCF (Discrete Correlation Function) 
method is used to analyze the correlation between the radio band and Gamma-ray and optical R band. 
It is found that there is a strong correlation between optical light variability and radio light variability, 


and the optical light variability is 66 + 40 days ahead of the radio variability. 
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